Contact graphs of disk packings as a model of spatial planar networks 
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Spatially constrained planar networks are frequently encountered in real-life systems. In this 
paper, based on a space-filling disk packing we propose a minimal model for spatial maximal planar 
networks, which is similar to but different from the model for Apollonian networks [J. S. Andrade, 
Jr. et al., Phys. Rev. Lett. 94, 018702 (2005)]. We present an exhaustive analysis of various 
properties of our model, and obtain the analytic solutions for most of the features, including degree 
distribution, clustering coefficient, average path length, and degree correlations. The model recovers 
some striking generic characteristics observed in most real networks. To address the robustness of the 
relevant network properties, we compare the structural features between the investigated network 
and the Apollonian networks. We show that topological properties of the two networks are encoded 
in the way of disk packing. We argue that spatial constrains of nodes are relevant to the structure 
of the networks. 

PACS numbers: 89.75.Hc, 89.75.Da, 05.10.-a 



I. INTRODUCTION 



The past decade has witnessed a great deal of activity devoted to complex networks by the scientific community, 
since many systems in the real world can be described and characterized by complex networks P, Q ■ Prompted by 
the computerization of data acquisition and the increased computing power of computers, researchers have done a lot 
. of empirical studies on diverse real networked systems, unveiling the presence of some generic properties of various 
natural and manmade networks: power-law degree distribution P{k) ~ k~ 7 with characteristic exponent 7 in the 
range between 2 and 3 [3[, small- world effect including large clustering coefficient and small average path length 
1^1 ' (APL) m, and degree correlations Q. These findings are important for our understanding of the real-life systems, 

since they strongly affect almost all aspects of various dynamical processes taking place on networks @, 0, [3 ■ 
O^l ' In order to reproduce or explain the above-mentioned striking common features of real-life systems, there has been 
^ \ a concerted effort in the last few years within the physical circle and elsewhere to develop network models to uncover 
and understand the complexity of real systems [l], 01 ■ In addition to the seminal Watts-Strogatz's (WS) small- world 
network model Q and Barabasi- Albert's (BA) scale- free network model Q, a huge variety of models and mechanisms 
. have been proposed to mimic real- world systems, including initial attractiveness @, aging and cost [13], fitness 
model [ll[ , duplication [l2| , weight or traffic driven evolution [13, 3, accelerating growth [TBI, [l6| , coevolution [I?} , 
visibility graph [IH, to name but a few. For reviews, see Refs. [H,0, sEl- Although significant progress has been made 
in the field of network modeling and has led to a significant improvement in our understanding of complex systems, it 
\ is still a fundamental task and of current interest to construct models mimicking real networks and reproducing their 
generic properties from different angles. 

Most previous network models concentrated on topological aspects and ignored the geographical effects. In these 
work, the node position has no special signification, which is reasonable for some networks that can be considered as 
lying in an abstract space, such as scientific collaboration network (l9j . metabolic network [2(jll and so forth. However, 
a plethora of real networks have well-defined node positions, including the Internet |2l|. power grid [22| . airline 
networks [23}, to name only a few. This class of networks is often referred to as geographical or spatial networks, 
which is a promising kind of networks, since its geography has a pivotal influence on the dynamics running on the 
networks, such as robustness (24|, cascading breakdown (25|, synchronization (26|, disease spreading J27[, among 
others. In addition to the geography, some real-world geographical networks are also planar. Typical examples are 
street networks [28j | , electronic circuit networks pfjj , ant-trail networks [30| , and neural networks [31 1 . Recently several 
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FIG. 1: The first two generations for the construction of the disk packing. 



models [13, [H, [H, HH, [35|, [3(| H3] have been presented that take the the geographical position into consideration. 
However, models for spatial planar networks are much less, with the exception of space-filling networks — Apollonian 
networks [38| that have received a great amount of attention [39|, HO, Hi], |42|, EH . 

In the paper, inspired by the disk packing and the construction method of Apollonian networks [38| . we propose a 
spatial planar network, where nodes (located at the center of disks) correspond to the disks and edges represent contact 
relation. The suggested network belongs to a deterministically growing class of self-similar graphs. On the basis of the 
recursive construction and self-similar structure of the network, we determine analytically many relevant topological 
properties, such as degree distribution, clustering coefficient, average path length, as well as degree correlations. The 
obtained results show that the network is simultaneously scale-free, small-world, and disassortativc, as observed in 
a variety of real networks. Except for the fact that all nodes of the network have separate spatial positions that 
encode the network structure, we also show that the network is not only planar, but also maximally planar. Note that 
although Apollonian networks are also derived from a recursive disk packing, a main subject of the current paper is 
to compare the topological features of the presented network and Apollonian networks, with an aim to address the 
effect of the way of disk packings (thus the spatial constrains of nodes) on the network structure, which is important 
to the theory of complex networks. 

II. CONSTRUCTION OF THE NETWORK 

Before defining the network, we first introduce a disk packing fiil . |4o| . which is a variation of the celebrated 
Apollonian packing [46| . The introduced disk packing, shown in figure [U is constructed as follows. We start with 
three mutually touching disks, the interstice of which is a curvilinear triangle, and we denote this initial configuration 
by generation t = 0. Then in the first generation t = 1, three smaller mutually touching disks are added to fill the 
interstice of the initial configuration, each of which touches two adjacent sides of the initial curvilinear triangle, giving 
rise to seven new smaller curvilinear triangles. For subsequent generations we indefinitely repeat the packing process 
for all the new curvilinear triangles. In the limit of infinite t generations, we obtain a disk packing. 

The above-mentioned disk packing can be used as a basis for a network, which is the research object of current 
paper. The translation from the disk packing to network generation is quite straightforward. Let the nodes (vertices) 
of the network correspond to the disks and make two nodes connected if the corresponding disks are in contact. 
Alternatively, one can also connect the centers of the touching disks by lines to obtain the network. Figure shows 
the network corresponding to the disk packing in figure [TJ 

In the construction process of the disk packing, for each interstice at arbitrary generation, once we add three disks 
to fill it, seven new interstices are created that will be filled in the next iteration. When building the network, it is 
equivalent to say that for each group of three new nodes added, seven new triangles are generated in the network, 
into each of which a group of three nodes will be inserted in the next iteration. According to this, we can introduce 
an iterative algorithm to create the network, denoted by F t after t generation evolutions. 



FIG. 2: Illustration of the network corresponding to the disk packing shown in figure [T] 




FIG. 3: Iterative construction method for the network. 



The iterative algorithm for the network is as follows: For t = 0, Fq consists of three nodes forming a triangle. Then, 
we add three nodes into the original triangle. These three new nodes are linked to each other shaping a new triangle, 
and both ends of each edge of the new triangle arc connected to a node of the original triangle. Thus we get Fi, see 
figure [3l For t > 1, F t is obtained from F t -\. For each of the existing triangles of F t -\ that has never generated 
nodes before, we call it an active triangle. We replace each of the existing active triangles of -Ft-i by the connected 
cluster on the right hand side of figure [3] to obtain Ft ■ The growing process is repeated until the network reaches a 
desired order (node number of network). Figure [5] shows the network growing process for the first two steps. 

Next we compute the order and size (number of all edges) of the network F t . Let L v (t), L e (t) and L&(t) be the 
number of nodes, edges and active triangles created at step t, respectively. By construction (see also figure |3J), each 
active triangle in Ft-x will be replaced by seven active triangles in Ft- Thus, it is not difficult to find the following 
relation: L&(t) = 7 L&(t — 1). Since La(0) = 1, we have L&{t) = 7*. 

Note that each active triangle in F t -\ will lead to an addition of three new nodes and nine new edges at step t, then 
one can easily obtain the following relations: L v (t) = 3L&(t — 1) = 3 x 7* _1 , and L e (t) = 9L&(t — 1) = 9 x 7* _1 for 
arbitrary t > 0. From these results, we can compute the order and size of the network. The total number of vertices 
N t and edges E t present at step t is 




(1) 



and 




(2) 



4 



respectively. So for large t, the average degree k t = -^p- is approximately 6, which shows the network is sparse as 
most real systems. 

From equations. {!]) and (|2|), we have E t = 3N t — 6. In addition, by the very construction of the network, it is 
obvious that arbitrary two edges in the network never cross each other. Thus our network is a maximal planar network 
(or graph) [47[ . 



III. STRUCTURAL PROPERTIES OF THE NETWORK 



In this section, we study the statistical properties of the network, in terms of degree distribution, clustering coeffi- 
cient, average path length, and degree correlations. 



A. Degree distribution 

When a new node i is added to the graph at step U (ti > 1), it has a degree of 4. Let Lj\(i,t) be the number of 
active triangles at step t that will create new nodes connected to node i at step t + 1. Then at step ti, L&(i, U) = 4. 
From the iterative generation process of the network, one can see that at any step each two new neighbors of i generate 
three new active triangles involving i, and one of its existing active triangle is deactivated simultaneously. We define 
ki(t) as the degree of node i at time t, then the relation between ki(t) and L&(i,t) satisfies: 

L A (i,t) = ki(t). (3) 

Now we compute L&(i,t). By construction, L&(i,t) = 3L&(i,t — 1). Considering the initial condition L^(i,ti) = 4, 
we can derive L^i^t) = 4 x 3 t_ti . Then at time t, the degree of vertex i becomes 

ki(t) = 4 x 3*-'\ (4) 

It should be mentioned that the initial three vertices created at step have a little different evolution process from 
other ones. We can easily obtain: L&(0, t) = 3* and ki(t) = 3* + 1. Thus, at step t, the degree of the initial three 
nodes is less than that of those three nodes born at step 1 but larger than that for those nodes emerging at other 
steps. 

Equation ffl shows that the degree spectrum of the network is discrete. It follows that the cumulative degree 
distribution |6[ is given by 

P m-V^l-I^ (51 

KnmW - 2^ N ~ 7t 5 , (5) 

T<ti 



which is valid for all ti > 2 

Substituting for ti (ti > 2) in this expression using ti = t j^y- gives 



When t is large enough, one can obtain 



7* x f^-l 1 " 7 / 1 ' 13 ) + 5 
P^(k) = U) ?t + 5 +5 - (6) 



, ^ - In 7/ In 3 



>(*) = I 4 ) • ( 7 ) 



So the degree distribution follows a power law form with the exponent 7=1 + |sZ , 



B. Clustering coefficient 

The clustering coefficient C l of node i is defined as the ratio between the number of edges that actually exist 
among the fej neighbors of node i and its maximum possible value, ki(ki — l)/2, i.e., C, = 2ei/[ki(ki — 1)]. The 
clustering coefficient of the whole network is the average of C i s over all nodes in the network. 

For our network, the analytical expression of clustering coefficient C(k) for a single node with degree k can be 
derived exactly. When a node enters the system, both ki and ej are 4. In the following iterations, each of its active 
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FIG. 4: Semilogarithmic plot of average clustering coefficient Ct versus network order Nt. 

triangles increases both ki and e, by 2 and 3, respectively. Thus, is equal to 4 + 1 (ki — 4) for all nodes at all steps. 
So one can see that there exists a one-to-one correspondence between the degree of a node and its clustering. For a 
node of degree k, we have 



C(k) 



2e _ 2 [4+ f(fc-4)] _ 4 1_ 

k(k- 1) _ k(k- 1) ~ k ~ k-l' 



(8) 



In the limit of large k, C(k) is inver sely proportional to degree k. The same scaling of C(k) ~ has also been 
observed in several real- life networks [481 ] . 

Using equation ((5J), we can obtain the clustering Ct of the network at step t: 



r=0 



Mr) /_4_ 
7V t VD r 



1 



D r - 1 



(9) 



where the sum runs over all the nodes and D r is the degree of the nodes created at step r, which is given by 
equation (j4]). In the infinite network order limit (Nt — > oo), equation ([9]) converges to a nonzero value C — 0.603, as 
shown in figure 2] Therefore, the average clustering coefficient of the network is very high. 



C. Average path length 

We represent all the shortest path lengths of network F t as a matrix in which the entry dij is the distance between 
node i and j that is the length of a shortest path joining i and j. A measure of the typical separation between two 
nodes in Ft is given by the average path length (also called average distance) dt defined as the mean of distances over 
all pairs of nodes: 

dt = Nt{Nt- l)/2 ' (10) 

where 

D t = Yl d » ( n ) 

denotes the sum of the distances between two nodes over all couples. Note that in equation (fTTj) . for a pair of nodes 
i and j (i ^ j), we only count d^ or dji, not both. 



FIG. 5: Schematic illustration of second construction means of the network. F t +i may be obtained by joining seven copies of 
Ft denoted as (77 = 1, 2, • • • , 7), which are connected to one another at the six edge nodes, i.e., A, B, C, X, Y, and Z. 



1. Recursive equation for total distances 

We continue by exhibiting the procedure of determining the total distance and present the recurrence formula, which 
allows us to obtain Dt+i of the t + 1 generation from Dt of the t generation. The network Ft under consideration has 
a self-similar structure that allows one to calculate D t analytically [i^, [5(| HH- As shown in figure EJ network F t+1 
may be constructed by joining at six edge nodes (i.e., A, B, C, X, Y, and Z) seven copies of F t that are labeled as 

According to the second construction method, the total distance D t+1 satisfies the recursion relation 

A+i = 7A + At-9, (12) 

(v) 

where At is the sum over all shortest path length whose endpoints are not in the same F t branch. The last term 
-9 on the right-hand side of equation (|12p compensates for the overcounting of certain paths: the shortest path cLax 
between A and X, with length 1, is included in both F^ and F^; similarly, the shortest paths cIay, dgx, dsz, 
dcY, dcz, dxY, dxz, and dyz are all computed twice. To determine D t , all that is left is to calculate A t . 

2. Definition of crossing distance 

In order to compute A t , we classify the nodes in F t +\ into two categories: the six edge nodes (such as A, B, C, 
X, Y, and Z in figure [5]) arc called connecting nodes, while the other nodes are named non-connecting nodes. Thus 
At, named the crossing distance, can be obtained by summing the following path length that are not included in 
the distance of node pairs in : length of the shortest paths between non-connecting nodes, length of the shortest 
paths between connecting and non-connecting nodes, and length of the shortest paths between connecting nodes (i.e., 
dAZ, d B Y, and d C z)- 

Denote A"' 13 as the sum of all shortest paths between non-connecting nodes, whose end-points are in F t {a) and 
.Fj . , respectively. That is to say, A"'^ rules out the paths with end-point at the connecting nodes belonging to F^ 
or F^ . For example, each path contributed to Aj' 2 does not end at node A, B, X or Y. On the other hand, let fij 



7 




FIG. 6: Illustration of the classification of interior nodes in F t (77 = 1, 2, • ■ • , 7), from which we can derive recursively the 
classification of interior nodes in network Ft+i. 



be the set of non-connecting nodes in F t +\, which belong to F t f . Then the total sum At is given by 

7 6 

A t = J2 E A t ,(3 + E d * + E *Bi + E d a + E dx ? 

P=a+la=l jen^ jen? 3 enJ jeoj 

??e{4,5,6,7} ??e{2,3,4,7} r f e{l,2,6,7} i)G{3,4,5} 



+ E dv J + E dz i + dAZ + dBY + d ° x ( 13 ) 

r ; e{l,5,6} )je{l,2,3} 

By symmetry, equation (|13[) can be simplified as 

At = 9 A, 1 ' 2 + 9 A 1 ' 3 + 3 A 1 ' 4 + 21 ^ d Aj + 6. (14) 

Having A t in terms of the quantities of A 1 ' 2 , A^' 3 , A 1,4 , and X^eo 4 ^Aj' ^ nc ucxt step is to explicitly determine these 
quantities. 



3. Classification of interior nodes 



To calculate the crossing distance A]' 



At 1 ' 3 , 



A 1 ' 4 , and £\ 



i t , aim ZsjEn* d'Aj, we classify interior nodes in network F t +\ into 
seven different parts according to their shortest path lengths to each of the three peripheral nodes (i.e. A, B, C). 
Notice that nodes A, B, C themselves are not partitioned into any of the seven parts represented as Pi, P2, P3, P4, 
P5, Pq, and P7, respectively. The classification of nodes is shown in figure For any interior node v, we denote the 
shortest path lengths from v to A, B, C as a, b, and c, respectively. By construction, a, 6, c can differ by at most 1 
since vertices A, B, C are adjacent. Then the classification function class(v) of node v is denned to be 



class(v) 



Pi 
Pi 
P>> 
Pi 
P5 
P<, 
Pi 



for a < b = c, 
for b < a = c, 
for c < a = 6, 
for a = c < 6, 
for a = b < c, 
for b — c < a, 
for a = b = c. 



(15) 
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It should be mentioned that the definition of node classification is recursive. For instance, class Pi and P4 in 



Ft belong to class P x in Ft+i, class P 3 and P 5 in P t w belong to class P 2 hi F t +i, class P 2 , Pq, and P7 in F t 
belong to class P5 in Pt+i. Since the three nodes A, B, and C are symmetrical, in the network we have the following 
equivalent relations from the viewpoint of class cardinality: classes Pi, P 2 , and P3 are equivalent to one another, and 
it is the same with classes P4, P5, and P§. We denote the number of nodes in network Ft that belong to class Pi 
as Nt,p 1 , the number of nodes in class Pi as Nt,p 2) and so on. By symmetry, we have Nt,p 1 = Nt y p 2 = Nt,p 3 and 
Nt p 4 = N t p 5 = N t p 6 . Therefore in the following computation we will only consider Ntp 1 , A^p 4 , and N t p 7 . It is 
easy to conclude that 



(i) 



N t = N t , Pl + N t ,p 2 + N t ,p 3 + N t ,p 4 + N t ,p 5 + N t ,p 6 + N t ,p 7 + 3 
= 3iV t: P 1 +3 7V t ,p 4 +^ f ,p 7 +3. 



(16) 



Considering the self-similar structure of the network, we can easily know that at time t + 1, the quantities Nt+i^, 
Nt+i.Pi, and N t +i t p 7 evolve according to the following recursive equations 



N t+ i. Pl = 3N t , Pl +4 N t , Pi + N ttPr , 
N t+ i.p 4 = 4N t , Pl + N t , Pi + N t , Pr + 1 , 
N t +i.P 7 = 6N t , Pi +N t:P7 , 



(17) 



where we have used the equivalent relations N tt p 1 = N t ,p 2 = N t ,p 3 and N t ,p i = N tt p 5 = N t ,p 6 . With the initial 
condition N 2 ,p 1 — 4, N 2 ,p 4 — 2, and N 2 ,p 7 = 6, we can solve the recursive equation (|T7|) to obtain 



Nt, Pl 
Nt, Pi 

N t ,Pr 



1 

124 



J_ 

02 



62 + 10 ■ 7* + 26 (-1 - a/2)' + 26 (-1 + v 7 ^)* + ll\/2 (-1 + n/2)' - llx/2 (-1 - V%) 
8 ■ 7* - 4 (-1 - y/Tf - 4 (-1 + V2) 1 + 15V2 (-1 + V2)* - 15^ (-1 - , 
- 4 • 7* - 33 (-1 - v/2)* - 33 (-1 + V2)* + 39\/2 (-1 - V2)' - 39V2 (-1 + V2) 



1 

124 



(18) 



62 



For a node i> in network Pt+i, we are also interested in the smallest value of the shortest path length from v to any 
of the three peripheral nodes A, B, and C . We denote the shortest distance as /.„, which can be defined to be 



f v = min(a, b, c). 



(19) 



Let dt ! p 1 denote the sum of /„ of all nodes belonging to class Pi in network F t . Analogously, we can also define the 
quantities dt,p 2 , d*,p 3 , ■ •-, dt,p 7 . Again by symmetry, we have dt,p 1 — dt,p 2 — dt,p 3 , dt,p 4 = <i*,p 5 = dt,p 6 , and dt,p r , 
dt,Pi, dt,p 7 can be written recursively as follows: 



= 3 d t ,p 1 + 4 d t ,p 4 + d t ,p 7 , 

= 4 d t , Pl + d t ,p 4 + d t .p 7 + 4 N t , Pl + 1 , 
= 6 (d t ,p 4 + N t ,p 4 ) + (d t ,p 7 + N t ,p 7 ) . 

Substituting equation (|18[) into equation (|20p , and considering the initial condition d 2 .p x = 4, c?2.p 4 
d 2 ,p 7 = 12, equation (|20|) is solved inductively 




(20) 



2, and 



H,P, 



H.Pi 



2 t,P 7 



5004888 



7 (-3 (-1 + v^* (20082 + 18235V2) + (-1 - \72)* (-60246 + 54705\/2) + 4 (29791 + 3327* 



5004888 



5004888 



+93 (l760 x 7* - 7 (-1 + x/2)' (1256 + 1149^/2) + (-1 - v^)* (-8792 + 8043^2)) t 
14 (-59582+ 13700 x 7* + (22941 - 56145^2) (-1 - V^)* + 3 (-1 + V2)* (7647+ 18715^2) 
+93 (l408 x 7* - 7 (-1 + V2)' (1042 + 107^2) + (-1 - x/2)* (-7294+ 749^2)) t , 
7 (9 (-1 - x/2)* (-16540 + 33393 V2) - 9 (-1 + V2)* (16540 + 33393\/2) + 8 (29791 + 7424 x 7*) 
+186 (704 x 7* - 21 (-1 - V2)' (-1149 + 628^2) + 21 (-1 + \/2)* (1149 + 628^)) t 



(21) 



4- Calculation of crossing distances 



Having obtained the quantities Nt p t and dt,p t (i = 1, 2, • • • , 7), we now begin to determine the crossing distance 
Aj' 2 , Aj' 3 , Aj' 4 , and X^jen 4 ^Aj expressed as a function of N t .p i and d t .p i . Here we only give the computation details 
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of A 4 ' 2 , while the computing processes of A* ' 3 , A 4 ' 4 , and X^jen 4 d A j are similar. For convenience of computation, we 
use r^' 1 to denote the set of interior nodes belonging to class Pi in ' . Then A 1,2 can be written as 

«erj>* , ie{i,2,3,4,5,6,7} 

The seven terms on the right-hand side of equation are represented consecutively as 5\ (i = 1, 2, • • • , 7). Next we 
will calculate the quantities 8\. By symmetry, 5 1 = 5 2 , = 5f . Therefore, we need only to compute 5\, Sf, Sf, 5% 
and Sj. Firstly, we evaluate 5j. By definition, 



51 J2 d - 



v^A,X,Y 

= Y ( d uA + d Av )+ Y (duA + d A Y + d Yv ) + 2J (duA + d A x + d Xv ) 

i€{i,4,6,7} ver 2 t - 3 vtr 2 t - 2 \jr 2 t - 5 

= N ttPl {'Sd t , Pl + 3d t ,P 4 + d tt p T + 2N t , Pl + N t ,Pt) + d t , Pl (3N ttPl + 3JV t ,p 4 + N t ,P 7 ). (23) 

Proceeding similarly, we obtain 

<5 3 = N t , Pl (3d t , Pl + 3d t ,p 4 + d t ,P 7 + 4N t , Pl + 3N t ,p 4 + JV t ,P 7 ) + d t , Pl (3N t , Pl + 3JV t ,p 4 + JV t ,P 7 ), (24) 

St = N tt p 2 (3dt, Pl + 3d t ,p 4 + d t , Pr + N t . Pl ) + d t ,p 2 (3Nt tPl + 3N tjP4 + N ti p 7 ), (25) 
<5 t 5 = N t ,p 2 (3d t , Pl + 3d t , Pi + d ti p 7 + 2JV t , Pl + N t .p 4 ) + d t , P2 (3N t , Pl + 3N t ^ + N t ,p 7 ), (26) 

and 

Sj = N tt P 3 (3dt, Pl + 3d t . Pi + d t ,P 7 + N t , Pl ) + d t ,P 7 (3N ttPl + 3N t .p 4 + N t>Pf ). (27) 
With the obtained results for 5 l t , we have 

A, 1 ' 2 = 2(3d t , Pl + 3d t , Pi + d tt P 7 )(3N tiPl + 3JV t ,p 4 + N t ,p,) + N t , Pl (3N t , Pl + 3N t ,P 4 + N t ,p 7 ) 

+2(N t , Pl + N t ,p 4 )(2N t , Pl + N t ,p 4 ) + N t , Pl (N t , P4 + N ttP7 ) + (JV t , Pl ) 2 . (28) 

Analogously, we find 

A 4 ' 3 = 2(3d t , Pl + 3d t ,P 4 + d t ,P 7 )(3N t , Pl + 3iVt,p 4 + JV t ,p 7 ) + 2(JV t , Pl ) 2 + 27V f , Pl (3A^ Pl + 3iV t ,p 4 + N t , P7 ) 

+N t , P4 (3N t , Pl + 3N t , Pi + N t ,p 7 ) + {N t , Pl + 2N t . Pi + N t ,p 7 ){2N t , Pl + N t , P4 ), (29) 

A, 1 ' 4 = 2(3d t , Pl + 3d t ,p 4 + d t ,p 7 )(3A^ Pl + 3iVt,p 4 + JV t ,p 7 ) + (3iV t , Pl + 3JV t ,p 4 + N t ^ 7 f + 3(JV t , Pl ) 2 , (30) 

and 

Y d Aj = (3dt.fi + 3dt.p 4 + d ti p 7 ) + (3JVt.fi + 3JV t , Pl + N t> p 7 ) + N tiPl . (31) 

Substituting equations (f2"8"|) , (j2"9")l . (j3"0|) , and (|3~lj) into equation HU), we the final expression for cross distances A t , 

1 



A 



3844 



5766 + 2728 x 7* + 30290 x 49* - 10974 (-1 + ^2 
9300\/2 (-1 +V2\ + 3114 (-1 + V2 
1017V2 (-1 + V2] 21 + 186 (-1 - V2V f -59 + 50^ 



9 (-1 - V2) 2t (-346 + 113V2) + 16368i x 49* . (32) 
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5. Exact result for average path length 



With the above-obtained results and recursion relations, we now readily calculate the sum of the shortest path 
lengths between all pairs of nodes. Inserting equation (j3"2")) into equation ([i"2")) and using the initial condition D 2 = 717, 
equation (|12p is solved inductively, 



D t = 



1 



161448 



201810 + 215270 x 7* + 11194 x 7 2t 



+ 72072 f-1 - y/2) - 57834\/2 (-1 - \/2) 



- 44037 (-1 - V2\* - 11340\/2(-l - y/fj* 
+ 72072 f-1 + V2j + 57834\/2 (-1 + \/2) 

- 44037 f-1 + V^) Vll340\/2(-l + \/2) ' 



16368* x 7* + 16368* x 7 2t 



(33) 



Substituting equation ()33|) into equation (|10|) yields the exactly analytic expression for average path length 

201810 



20181(15 + 8 x 7' + 7 2t ) 
215270 x 7* + 11194 x 7 2t 

72072 (-I-V2) - 57834 V2 f-1 - 



44037 (-1 - V2) 2t - 11340^ f-1 - v^V 



72072 (-I + V2) 1 + 57834 V2 f-1 + v^V 



44037 (-1 + V2) 2 ' + 11340\/2 (-1 + V2 



2/ 



16368* x 7* + 16368* x 7 2t 



(34) 



In the large * limit, d t ~ *, while the network order N t ~ 7* which is obvious from equation (fT]). Thus, the 
average path length grows logarithmically with increasing order of the network. We have checked our analytic result 
provided by equation (|34[) against numerical calculations for different network order up to < = 8 which corresponds 
to N% = 1007 772. In all the cases we obtain a complete agreement between our theoretical formula and the results 
of numerical investigation, see figure [7J 

Recently, it has been suggested that for random uncorrelated scale-free networks (SFNs) with degree exponent 
7 < 3 and network order N, their average distance d(N) behaves as a double logarithmic scaling with N: d(N) ~ 
In In N [52, 53]. However, for the deterministic network considered here, in despite of the fact that its degree exponent 
7 = 1 + < 3, its average path length scales as a logarithmic scaling with network order, showing a obvious 
difference from that of the stochastic scale-free counterparts. The logarithmic scaling of dt with Nt for our network 
as well as the Apollonian networks [4l[ shows that previous relation between APL and the network order obtained 
for uncorrelated SFNs [13, [53[ is not valid for disassortative SFNs [52|, [53[, at least for some spatial networks, e.g., 
the Apollonian networks and the network considered here. This leads us to the conclusion that degree exponent itself 
does not suffice to characterize the APL of SFNs. 



D. Degree correlations 

An interesting quantity related to degree correlations is the average degree of the nearest neighbors for nodes with 
degree k, denoted as k nn (k), which is a function of node degree k [5J, |55j. When k nn (k) increases with k, it means 
that nodes have a tendency to connect to nodes with a similar or larger degree. In this case the network is defined 
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N 
t 

FIG. 7: Average distance d t versus network order N t on a semilogarithmic scale. The solid line serves as guide to the eye. 



as assortative In contrast, if k nn (k) is decreasing with k, which implies that nodes of large degree are likely 
to have near neighbors with small degree, then the network is said to be disassortative. If correlations are absent, 
kuri(k) = const. 

We can exactly calculate k nn (k) for network F t using equations ([3]) and ^ to work out how many links are made 
at a particular step to nodes with a particular degree. By construction, we have the following expression (56l . \Et\ 



k n n{k) 



L v (U)k(U,t) 



t'=t 



*i=t*-i 

J2 2^(^(4^-1^(4*) 



t'=o 



2L v (U)L A (U,t' t -l)k(t' l ,t) 



t'=ti+i 



+ 2 



(35) 



for k = 4 x 3*~ 4i (ti > 1), where k(ti,t) is the degree of a node i at time t that was born at step ti. Here the first 
sum on the right-hand side accounts for the links made to nodes with larger degree (i.e. t[ < U) when the node was 
generated at ti. The second sum describes the links made to the current smallest degree nodes at each step t' i > ti. 
The last term 2 accounts for the two links connected to two simultaneously emerging nodes. After some algebraic 
manipulations, we can rewrite equation (|35p in term of k to obtain 



knn(k) 



(3 



2t+l 



For k = 3* + 1 (U = 0), we have 



k n n{k) — 



' 4 x2-ln 7/ ln3 



2x7 



t-1 



8M|) 
31n3 



10. 



(36) 



1 



k(t h t) 

8t x 3* 
3 <+1 + 3 



2L A {t i ,t' i -l)k{t' i ,t) 



t'.=ti+i 
+ 2. 



(37) 



Therefore, for large t and k, k nn (k) is approximately a power law function of k as fc nn (fc) ~ k~ u with uj = 2 — ~ 
0.229, which shows that the network is disassortative. Note that k nn (k) of the Internet exhibits a similar power-law 
scaling with exponent u> = 0.5 [54| . 



IV. CONCLUSION 



In summary, motivated by the disk packing and Apollonian networks, we have presented a model for spatial planar 
networks introducing the influence of geography encoded in the disk packing. According to the construction, we have 
studied analytically the main structural features of the network. We have shown that the network has a power law 
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distribution with exponent 7=1 + it has a large clustering coefficient 0.603, its APL scales logarithmically with 
the number of network nodes, and it is disassortative with the average degree of the nearest neighbors for nodes 
having degree k being roughly a power-law function of k with exponent -0.229. 

Note that although both the network considered here and the Apollonian network [38|, [56| are translated from disk 
packings, and both networks have qualitatively similar topologies, their structural characteristics are quantitatively 
different. For example, the exponent of degree distribution is 1 + for our network, while for Apollonian network it 
is 1 + ■ the average clustering coefficients for our network and Apollonian network are 0.603 and 0.828, respectively. 
In addition, the average path length [4l| and the degree correlations [56| for both networks are also of quantitative 
difference. These disparities of the two networks show that the ways of disk packing lead to different spatial constrains 
of network nodes, which in turn have a significant impact on network properties and thus dynamics running on 
networks, such as cascading failing 25 1, random walks (58j . and so on. Thus, we can conclude that for spatial 
networks the positions, where nodes are geographically located, matter greatly and should be incorporated when 
modeling such networks. Ignoring the geography will lead to miss some important attributes and properties of the 
systems. 
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